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We present a measurement of the decay B~ r~V T using a data sample containing 657 x 10 6 
BB pairs collected at the T(4S) resonance with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider. A sample of B + B~ pairs are tagged by reconstructing one B + meson decaying 
semileptonically. We detect the B~ — > t~v t candidate in the recoil. We obtain a signal with a 
significance of 3.6 standard deviations including systematic uncertainties, and measure the branching 
fraction to be B(B~ — > t~V t ) — [l^lg'^stat) lo'g^syst)] x 10 -4 . This result confirms the evidence 
for B~ — > t~V t obtained in a previous Belle measurement that used a hadronic B tagging method. 

PACS numbers: 13.20.He, 14.40.Nd 



The purely leptonic decay B~ — > r~7/ r [l| is of par- 
ticular interest since it provides a unique opportunity to 
test the Standard Model (SM) and search for new physics 
beyond the SM. In the SM, the branching fraction of the 
decay B~ — > t~v t is given by 

nCR-^n--r,\ G F m Bm 2 T f m£\ 2| , 2 

B(B ->• t v T ) = 1 T f B \Vub\ T B , 

8tt V m B J 

(1) 

where Gf is the Fermi coupling constant, m T and tub are 
the r lepton and B~ meson masses, t b is the B~ life- 
time, \V u b\ is the magnitude of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element Q, and /g is the B 
meson decay constant. Dependence on the lepton mass 
arises from helicity conservation, which suppresses the 
muon and electron channels. A recent SM estimation of 
the branching fraction [3] is (0.76±qH) x 1(T 4 . In the 
absence of new physics, measurement of the B~ — > t~V t 
decay can provide a direct experimental determination 
of Jb, which can be compared to lattice QCD calcula- 
tions 4]. Physics beyond the SM. however, could signifi- 
cantly suppress or enhance B(B~ — > t~V t ) via exchange 
of a new charged particle such as a charged Higgs boson 
from supersymmetry or two-Higgs doublet models Q . 

Belle previously reported Q the first evidence of B~ — > 
t~v t decay with a significance of 3.5 standard devia- 
tions (er), and measured the branching fraction to be 
B(B- -> r-V T ) = (1.79t^(stat)±°f ( syst )) x lO" 4 , 
using a hadronic reconstruction tagging method. The 
BaBar Collaboration reported a search for B~ — > t~v t 
decay with hadronic tagging [8(] using 383 x 10 6 BB pairs 
and with semileptonic tagging Q using 459 x 10 6 BB 
pairs. Combining the two measurements, they obtained 



a 2.8cr excess and a branching fraction B(B~ — s- t~77 t ) = 
(1.7±0.6) x 1CP 4 . These experimental results are slightly 
larger than the SM estimation in Ref. [3|, though the 
statistical precision is not sufficient to demonstrate the 
existence of physics beyond the SM. To better estab- 
lish this decay mode and determine the branching frac- 
tion with greater precision, we present a measurement 
of B~ t~v t from Belle using a semileptonic tagging 
method. 

We use a 605 fb _1 data sample containing 657 x 10 6 
BB pairs collected with the Belle detector at the KEKB 
asymmetric-energy e + e~ (3.5 on 8 GeV) collider [l(| op- 
erating at the T(45) resonance (y/s = 10.58 GeV). We 
also use a data sample of 68 fb~ 4 taken at a center of 
mass energy 60 MeV below the nominal T(4S) mass 
(off-resonance) for background studies. The Belle detec- 
tor is a large-solid- angle magnetic spectrometer that 
consists of a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel thresh- 
old Cherenkov counters (ACC), a barrel-like arrangement 
of time-of-flight scintillation counters (TOF), and an 
electromagnetic calorimeter (ECL) comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux- return 
located outside of the coil is instrumented to detect K\ 
mesons and to identify muons (KLM). Two inner detector 
configurations were used. A 2.0 cm beampipe and a 3- 
layer silicon vertex detector were used for the first sample 
of 152 xlO 6 BB pairs, while a 1.5 cm beampipe, a 4-layer 
silicon detector and a small-cell inner drift chamber were 
used to record the remaining 505 xl0 6 i?B pairs [Hj]. 

We use a detailed Monte Carlo (MC) simulation based 
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on GEANT [13j to determine the signal selection effi- 
ciency and study the background. In order to reproduce 
the effects of beam background, data taken with random 
triggers for each run period are overlaid on simulated 
events. The B~ — > t~V t signal decay is generated by 



the EvtGen package [14[. Radiative effects are modeled 
using the PHOTOS code 15j. To model the background 



from 



— > BB and continuum qq (q = u,d,s,c) pro- 



duction processes, we use large MC samples of BB me- 
son pair decays to charm and continuum qq processes 
corresponding to about ten times and six times the data 
sample, respectively. We also use MC samples of rare 
B decay processes such as charmless hadronic, radiative, 
electroweak decays and b — > u semileptonic decays. The 
contamination from other low multiplicity backgrounds 
such as e + e~ — > t + t~ and two-photon processes is also 
studied using dedicated MC samples. 

The B~ — » t~77 t candidate decays are selected us- 
ing the feature that at the T(AS) resonance B meson 
pairs are produced with no additional particles. We first 
reconstruct one of the B mesons decaying semileptoni- 
cally (referred to hereafter as Stag) and then compare the 
properties of the remaining particle(s) in the event (B S j g ) 
to those expected for signal and background. In order to 
avoid experimental bias, the signal region in data is not 
examined until the event selection criteria are finalized. 

Charged particles are selected from well measured 
tracks (reconstructed with the CDC and SVD) originat- 
ing from the interaction point. Electron candidates are 
identified based on a likelihood calculated using the fol- 
lowing information: dE/dx measured in the CDC, the 
response of the ACC, the ECL shower shape and the ra- 
tio of the ECL energy deposited to the track momentum. 
Muon candidates are selected using KLM hits associated 
to a charged track. Both muons and electrons are se- 
lected with efficiency greater than 90% in the momen- 
tum region above 1.2 GeV/c, and misidentification rates 
of less than 0.2% (1.5%) for electrons (muons). After 
selecting leptons, we distinguish charged kaons from pi- 
ons based on a kaon likelihood derived from the TOF, 
ACC, and dE/dx measurements in the CDC. The typi- 
cal kaon identification efficiency is greater than 85% and 
the probability of misidentifying pions as kaons is about 
8%. Photons are identified as isolated ECL clusters that 
are not matched to any charged track. Neutral tt° can- 
didates are selected from pairs of photons with invariant 
mass between 0.118 and 0.150 GeV/c 2 . The energy of the 
photon candidates must exceed: 50 MeV for the barrel, 
100 MeV for the forward endcap and 150 MeV for the 

backward endcap, except for low momentum n° candi- 
— *o — o n 

dates from D — > D tt decay for which we require the 
photon energy to be greater than 30 MeV. 

We reconstruct the i? tag in B + — » D £ + vi and B + — » 
D°£ + ise decays, where £ is electron (e) or muon (/i). D 
mesons are reconstructed in the K + ir~ , K + tt~it and 



K + n~TT + n~ modes. For B s ; g , we use r~ decays to only 
one charged particle and neutrinos i.e. r~ —¥ l~vgv T and 

T~ — > 1T~V r . 

We require the invariant mass of D candidates 
to be in the range [1.851 GeV/c 2 , 1.879 GeV/c 2 ] 



for D 



K + ir and K + tt 



decays, and 



[1.829 GeV/c 2 , 1.901 GeV/c 2 ] for D 
— *o 

decay. D candidates are selected by combin- 
ing the D candidates with low momentum it can- 
didates or photons. For D candidates, we re- 
quire the mass difference AM = M D ,o — M D a to 
be in the range [0.1389 GeV/c 2 , 0.1455 GeV/c 2 ] and 



[0.123 GeV/c 2 , 0.165 GeV/c 2 ] for D 



*° n° o a 
D it and 



D — > D 7 decays, respectively. These regions corre- 
spond to three standard deviations in the correspond- 

— *o 

ing resolutions. To suppress D 's from continuum back- 

— *o 

ground processes, the momentum of D candidates cal- 
culated in the T(4£>) center-of-mass system (cms) is re- 
quired to be less than 2.5 GeV/c. 

We select signal candidates from events with one D or 
D and one £ + to form -B tag , and one £ or it candidate 
for -B s ig- We require that no other charged particle or it 
remain in the event after removing the particles from 
the Stag and B s - lg candidates. The -Btag candidates are 
selected using the lepton momentum in the cms frame, 
P £ cms , and the cosine of the angle between the direction of 
the Stag momentum and the direction of the momentum 

(*)0 

sum of the D and the lepton, cos^ £>(*)£, measured 
in the cms frame. This angle is calculated using 



cos( 



B,D(')e 



O 77^cms 7?cms 



2 _ Ji,f2 

B Ivl D(')i 



o perns pems 



(2) 



where -E£™ s m is the beam energy, Pg ms is the momentum 
of B meson calculated with PJ, ms = V^bcam) 2 ~ m %> 
P^ )l and M D {,) 1 are the energy sum, momen- 

— (*)o 

turn sum and invariant mass, respectively, of the D 
and lepton system. All parameters are calculated in the 
cms. Properly reconstructed -Btag candidates are popu- 
lated within the physical range [—1,1], while combina- 
torial backgrounds can take unphysical values. For the 
signal side, the £~ or ir~ candidate from the r decay is 
selected using the momentum in the cms, denoted -P s C jg ls . 
The signal yield is obtained by fitting the distribution of 
the remaining energy in the ECL, denoted -Becl, which is 
the sum of the energies of ECL clusters that are not asso- 
ciated with particles from the -B tag and -B S i g candidates; 
here the -Eecl clusters satisfy the same minimum energy 
requirements as photon candidates. For signal events, 
.Eecl must be either zero or a small value arising from 
splitoff showers created by -Bt ag and B s i g particles and 
residual beam background hits. Therefore, signal events 
peak at low Becl- On the other hand, background events 
are distributed toward higher Becl due to the contribu- 
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tion from additional particles. The selection criteria for 
P^ cms , cos9 BD (,) t and P s c ™ s are optimized for each of the 
r decay modes, because the background levels and the 
background components are mode-dependent. The op- 
timization is done so that the figure of merit s/y/s~+~n 
is maximized, where s and n are the number of signal 
and background events expected in the signal-enhanced 
region Pecl < 0.2 GeV, calculated assuming a signal 
branching fraction of 1.79 x 10 -4 . For leptonic r decays, 
the dominant background is from PP events tagged by 
a semileptonic decay with a correctly reconstructed com- 

bination of a D and a t + . For these decays loose selec- 
tion criteria are chosen to maintain high signal efficiency: 
0.5 GeV/c < P/ ms < 2.5 GeV/c, -2.1 < cos8 b ,dh < 1-3 
for the D mode or —2.6 < cos9b,di < 1.2 for the D 
mode, and 0.3 GeV/c < P s °i™ s . For the hadronic r decay 
mode, there is more background from e + e~ — s- qq contin- 
uum and combinatorial D^°£ background. Tighter cri- 
teria are used to reduce such backgrounds: 1.0 GeV/c < 
P/ ms < 2.2 GeV/c, -1.1 < cos9 BMi ,)o e < 1.1, and 
1.0 GeV/c < P s c ™ s < 2.4 GeV/c. The upper bound 
on P s c ;™ s is introduced to reject two-body B decays. In 
addition, we suppress continuum background by requir- 
ing the cosine of the angle between the signal side pion 
track and the thrust axis of the B tag , cos6*thr, to be 
less than 0.9. We select candidate events in the range 
Pecl < 1-2 GeV for further analysis. The number of 
candidate events are 2481 for r~ — > e~v e v T: 2011 for 
t~ — > pTv il y T and 1018 for t~ — > ~k~v t decays. FigureQ] 
shows the cos 9 B £>(*)^ distribution for the signal candi- 
date events including both leptonic and hadronic r decay 



modes with all selection criteria other than cos( 
applied. The excess over the MC expectation for events 
without B — > D^lv decays indicates that the final sam- 
ple contains candidate events with the correct combina- 

tion of a D and a £ + forming a B tag . In the remaining 
candidates, according to a MC study, 4.6%, 13.4% and 
12.0% are events without a P tag from B + B~, B°B° and 
non-Pi? processes, respectively. 

The number of signal events is extracted from an ex- 
tended maximum likelihood fit to the Pecl distribution 
of the candidate events. Probability density functions 
(PDFs) for each r decay mode are constructed from the 
MC simulation. We use Pecl histograms obtained from 
MC samples for each of the signal and the background 
components. The PDFs are combined into a likelihood 
function, 



AH 



(3) 



= 1 3 



where j is an index for the signal and background con- 
tributions, rij and fj are the yield and the PDF, respec- 
tively, of the jth component, P; is the Pecl value in the 
ith event, and N is the total number of events in the 




B-D(*)l 



FIG. 1: cos 6 g distribution for candidate events with 

Pecl < 1.2 GeV selected with all Ptag and P s i g requirements 
except for those on cos6 B D {,) 1 . Leptonic and hadronic r 
decay modes are combined. The points with error bars are 
data. The dot-dashed, dotted and dashed histograms are the 
MC expectation for events without B + — > i + ve decays 
for P+B", sum of P+B~ and P°P°, and sum of BP and 
non-BP events, respectively. 



data. The dominant background components are from 
PP decays to a final state with charm and continuum 
processes. The small background from rare charmless B 
decays and other low multiplicity processes such as r pair 
and two-photon processes is also included in the fit. In 
the final sample with Pecl < 1-2 GeV, the fractions of 
the background from rare charmless P decays and low 
multiplicity non-P processes are estimated from MC to 
be 8% and 3% for leptonic r decays and 11% and 8% for 
hadronic r decay, respectively. 

The Pecl estimation in MC is validated using vari- 
ous control samples. The MC distributions of not only 
Pecl but also P £ cms , cos9 B D {, )ll P s c ;™ s and cos# t hr are 
compared to those of the control samples to confirm that 
MC describes the background composition properly. The 
off-resonance data is used to calibrate the MC simulation 
of the continuum background. We find that our MC un- 
derestimates the overall normalization of the continuum 
background though the predicted shapes are consistent 
with data within statistical errors. We obtain the cor- 
rection factor for the overall normalization of the con- 
tinuum MC to be 1.43 ±0.11 by comparing the number 
of remaining events in off-resonance data with the MC 

expectation. The sidebands in cos B D (,)£, D mass, the 
— *o — o 

mass difference between D and D , and Pecl are used 
as control samples to check the overall background de- 
scription including the PP contribution. The distribu- 
tions in these variables obtained from MC with the con- 
tinuum normalization correction applied are found to be 
consistent with the corresponding distributions in data. 
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(GeV) 



FIG. 2: -Becl distribution for double semileptonic tagged 
events. The points with error bars are data and the solid 
histogram is the MC expectation scaled to the luminosity of 
the data. The dashed histogram is the MC expectation mul- 
tiplied by the correction function described in the text. The 
inset shows the ratio of data to the MC expectation and the 
correction function. 



The agreement between MC and data is also confirmed 
in 73° tagged events where the B tag is reconstructed in 
E° —> D*~l + v decays. The contributions to the Eecl 
distribution are not only from beam background but also 
include splitoff showers originating from B tag and B slg 
decay products. The relative fractions of these sources 
are 21%, 53% and 26%, respectively, in the signal MC 
sample. To take into account the possible difference 
between MC and data descriptions of splitoff showers, 
the signal Eecl shape is calibrated using double tagged 
events, in which the -B tag is reconstructed in a semilep- 
tonic decay as described above and B s i g is reconstructed 
in the decay chain, B~ -> D*°£-V (D*° -> E>°7r°), fol- 
lowed by D° —> K~tt + . Figure [2] shows the Eecl dis- 
tribution in the double tagged sample for data and for 
the MC simulation scaled to the same luminosity. The 
background in this control sample is negligibly small. 
We find the -Eecl distribution of data tends to have a 
slightly smaller width than MC. The difference between 
the data and MC is parameterized as a first-order poly- 
nomial function of Eecl obtained by fitting the ratio of 
data to MC for the Eecl histograms of the double tagged 
sample. The ratio and the fit result are also shown in 
Fig. [2] The Eecl histogram obtained from the signal 
MC sample is multiplied by this correction function. 

In the final fit, four parameters are allowed to vary: 
the total signal yield and the sum of BB and contin- 
uum backgrounds for each r decay mode. The ratio 
of the BB to the continuum background is fixed to the 



value obtained from MC with the normalization correc- 
tion applied. Other background contributions are fixed 
to the MC expectation. We combine r decay modes by 
constraining the ratios of the signal yields to the ratio 
of the reconstruction efficiencies obtained from MC in- 
cluding the branching fractions of r decays [16|. Fig- 
ure [3] shows the Eecl distribution overlaid with the fit 
results. The Eecl distribution for each r decay mode 
is also shown. We see a clear excess of signal events 
in the region near zero and obtain a signal yield of 
n s = 143+|f. The branching fraction is calculated as 
B = n s /(2eN B + B -), where e is the reconstruction effi- 
ciency including the branching fraction of the r decay 
mode and N B + B - is the number of T(45) — > B + B~ 
events, assuming N B + B - = N ga -gO. Table Q] lists the 
signal yields and the branching fractions obtained from 
separate fits to each r decay mode and the fit with all 
three modes combined. The results of the individual fits 
are consistent within statistics. The \ 2 of the three re- 
sults is 2.43 for two degrees of freedom, corresponding to 
a 30% confidence level. 




FIG. 3: -Eecl distribution of semileptonic tagged events 
with the fit result for (a) all r decay modes combined, (b) 
r~ — > e~V<aj T , (c) — > fi^V^Ur and (d) t~ — > n~i/ T . The 
points with error bars are data. The hatched histogram and 
solid open histogram are the background and the signal con- 
tributions, respectively. 



Systematic errors for the measured branching fraction 
are associated with the uncertainties in the signal yield, 
efficiencies and the number of B + B~ pairs. Unless ex- 
plicitly stated otherwise, the systematic errors for each 
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TABLE I: Results of the fit for signal yields and branch- 
ing fractions, e is the reconstruction efficiency including the 
branching fraction of the r decay mode. The first error in the 
branching fraction is statistical and the second is systematic. 



Decay Mode 


Signal Yield 


e, 10" 4 


B, 1(T 4 


r~ —¥ e~V e v T 

T~ — > yTV^VT 
T~ —¥ TT~ V T 


70+23 

' °-22 
12 +18 
lz -17 
cc + 21 
dJ -20 


5.9 
3.7 
4.7 


1 Qn +0. 59+0.33 
± - yu -0.57-0.35 
n cn+0.76+0.18 
U.OU_ .72-0.21 
1 on+0-69+0.36 
i - ou -0.66-0.37 


Combined 


143li 


14.3 


1 r ,1+0.38 +0.29 
J -*"*— 0.37— 0.31 



source are obtained by varying the corresponding param- 
eters individually by their uncertainties, repeating the fit 
procedure and adding differences from the nominal re- 
sult in quadrature. The systematic errors for the signal 
yield arise from the uncertainties in the PDF shapes for 
the signal and for the background. The uncertainty in the 
signal shape correction function is estimated by changing 
the parameters of the correction function by their errors 
and replacing the function with a second-order polyno- 
mial (ij '4%)- The systematic error from MC statistics is 
evaluated by varying the content of each bin in the sig- 
nal -Eecl PDF histograms by its statistical uncertainty 
(±0.9%). The main contributions to the systematic er- 
rors for the background PDF shapes are statistical er- 
rors in the MC histograms which is estimated 
in the same way as the signal PDF MC statistical un- 
certainty. Other large sources are the uncertainties in 
the background composition. The errors due to the un- 
certainties in the branching fractions of B decay modes 
that peak near zero -Eecl such as B~ — > D°£~Vi with 
D° -+ K° L Kl, Klir" and K~l + v t , and B° -+ D + l~V t 
with D + — > K^(. + h>i are estimated by changing the 
branching fractions in MC by their errors [l6| (lg'f%). 
For branching fractions of D decays with a K\, we use 
the values for the corresponding D decays with Kg's. 
Uncertainties in the background from the possible contri- 
bution of rare charmless B decays such as B~ — > tt°£~V^, 
B~ — > K~vV and £~77^7, and from t + t~ pair and two 
photon processes are evaluated by changing the fractions 
obtained from MC by their experimental errors 16] if 
available, or by ±50% otherwise The system- 
atic error due to the uncertainty in the normalization 
correction factor for the continuum MC is The 
systematic error associated with the reconstruction ef- 
ficiency of the tag-side B is evaluated by comparing of 
the B(B~ — » D*°£~T7i) branching fraction measured with 
the double tagged sample in data to the world average 
value [lj|. We obtain the ratio to be 0.907 ± 0.044 and 
take the difference from unity plus one a as the system- 
atic error (13.7%). The systematic errors in the signal- 
side efficiencies arise from the uncertainty in tracking ef- 
ficiency (1.0%), particle identification efficiency (1.3%), 
branching fractions of r decays (0.4%), and MC statistics 



(0.8%). The systematic error due to the uncertainty in 
Nb+b- is 1-4%. The total fractional systematic uncer- 
tainty is 1 20 %i ano - the branching fraction is 

B(B- -> t-77 t ) = (1.54t°-3 7 (stat)±°-^(syst)) x 10~ 4 . 

(4) 

The significance of the observed signal is evaluated by 
E = y— 2 ln(£o/£max) where £ max and £0 denote the 
maximum likelihood value and likelihood value obtained 
assuming zero signal events, respectively. The systematic 
uncertainty is convolved in the likelihood with a Gaussian 
distribution having a width corresponding to the system- 
atic error of the signal yield. We find the significance of 
the signal yield to be 3.6a. 

In summary, we have measured the decay B~ — > t~77 t 
with BB events tagged by semileptonic B decays using 
a data sample containing 657 x 10 6 BB pairs collected 
at the T(45 l ) resonance with the Belle detector at the 
KEKB asymmetric-energy e + e~ collider. We measure 
the branching fraction to be (1.54lo;|?(stat) Ig'^syst)) x 
10~ 4 , with a significance of 3.6 standard deviations in- 
cluding systematics. This result is consistent with the 
previous Belle measurement using BB events tagged by 
hadronic B decays and is consistent with the results re- 
ported by the BaBar collaboration. Using the measured 
branching fraction and known values of Gf, tub, Tn T and 
t~b [IH , the product of the B meson decay constant fs 
and the magnitude of the Cabibbo-Kobayashi-Maskawa 
matrix element \V u b\ is determined to be 

f B \V ub \ = (9.311;? ± 0-9) x 10~ 4 GeV. (5) 

Using \ V ub \ = (3.89 ±0.44) x 10~ 3 in Ref. 16], f B is cal- 
culated to be 0.24 ± 0.05 GeV. The measured branching 
fraction is consistent within errors with the SM expecta- 
tion from other experimental constraints Q. The result 
can be used to extract constraints on new physics models. 

We thank the KEKB group for excellent operation 
of the accelerator, the KEK cryogenics group for effi- 
cient solenoid operations, and the KEK computer group 
and the Nil for valuable computing and SINET3 net- 
work support. We acknowledge support from MEXT, 
JSPS and Nagoya's TLPRC (Japan); ARC and DIISR 
(Australia); NSFC (China); MSMT (Czechia); DST (In- 
dia); MEST, NRF, NSDC of KISTI, and WCU (Korea); 
MNiSW (Poland); MES and RFAAE (Russia); ARRS 
(Slovenia); SNSF (Switzerland); NSC and MOE (Tai- 
wan); and DOE (USA). 



[1] The charge-conjugate decays are implied throughout this 

paper unless otherwise stated. 
[2] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 

652 (1973); N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). 



7 



[3] J. Charles et al. (CKMfitter Group) Eur. Phys. J. C 
41, 1 (2005); Preliminary results as of ICHEP 2010, 
|http://ckmfitter in2p3.fr/plotsJCHEP10 

[4] E. Gamiz, C. T. H. Davies, G. P. Lepage, J. Shigemitsu 
and M. Wingate (HPQCD Collaboration) Phys. Rev. D 
80, 014503 (2009); S. Aoki et al. (JLQCD Collaboration) 
Phys. Rev. Lett. 91, 212001 (2003); C. Bernard et al. 
(MILC Collaboration) Phys. Rev. D 66, 094501 (2002); 
A. Ali Khan et al. (CP-PACS Collaboration) Phys. Rev. 
D 64, 054504 (2001). 

[5] W. S. Hou, Phys. Rev. D 48, 2342 (1993). 

[6] S. Baek and Y. G. Kim, Phys. Rev. D 60, 077701 (1999). 

[7] K. Ikado et al. (Belle Collaboration) Phys. Rev. Lett. 97, 
251802 (2006). 

[8] B. Aubert et al. (BaBar Collaboration) Phys. Rev. D 77, 

011107(R) (2008). 
[9] B. Aubert et al. (BaBar Collaboration) Phys. Rev. D 81, 



051101(R) (2010). 

[10] S. Kurokawa and E. Kikutani, Nucl. Instrum. Methods 
Phys. Res., Sect. A 499, 1 (2003), and other papers in- 
cluded in this volume. 

[11] A. Abashian et al. (Belle Collaboration), Nucl. Instrum. 
Methods Phys. Res., Sect. A 479, 117 (2002). 

[12] Z. Natkaniec et al. (Belle SVD2 Group), Nucl. Instrum. 
Methods Phys. Res., Sect. A 560, 1 (2006). 

[13] R. Brun et al, GEANT3.21, CERN Report DD/EE/84-1 
(1984). 

[14] D. J. Lange, Nucl. Instrum. Methods Phys. Res., Sect. A 

462, 152 (2001). 
[15] E. Barbelio and Z. Was, Comput. Phys. Commun. 79, 

291 (1994). 

[16] K. Nakamura et al, J. Phys. G 37, 075021 (2010). 



